A0-A0a9  M6  WISCONSIN  UNIV-MAOISON  MATHEMATICS  RESEARCH  CENTER  F/6  12/1 

THEORETICAL  AND  PRACTICAL  ASPECTS  OF  SINBULARITT  AND  EieCNHOOE  — <TC(U> 
MAR  60  A  6  RAMM  DAA6-29-75-C-002R 

UNCLASSIFIED  MRC-TSR-20S1  NL 


^0  A089588 


Mathematics  Research  Center 
University  of  Wisconsin-Madison 
610  Walnut  Street 
Madison,  Wisconsin  53706 


March  1980 


(Received  February  25,  1980) 


U.  S.  Army  Research  Office 
p.  0.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 

I _ 


Approved  for  public  reloase 
Distribution  unlimited 


Air  Force  Office  of  Scientific  Research 
Bolling  AFB 
Washington,  EJC  20332 

80  8  11  066 


SIGNIFICANCE  AND  EXPT.ANATION 

In  ret.^=n^  '.’ecirs  a  vast  amount  of  literaturp  has  bepn  written  on 
sinqularitv  and  eiqonmodp  expansion  methods  (SEM  and  EEM) .  The  main  practical 
[-'■urriosr  is  to  ;  rodnce  a  theory  which  allows  us  to  identify  a  flyinq  (or 
stationary  v.;r n  t)  from  the  observed  transient  field  scattered  by  the 
tarqet.  ini.-;  ;  rdiloni  is  discussed  in  the  present  report.  We  describe  11 
tile  l^asi'  St  irtinq  coinrs  of  the  cnqinecrinq  ainiroach  to  the  problem,  and 
tile  extent  to  wiii,' h  they  are  consistent  f rtim  the  mathematical  point  of 
view;  2)  tlio  material  which  has  been  riqcronsiv  established  in  the  field 
by  trie  writer  and  other  authors;  1)  the  imi>ort.,mt  I'oints  in  jitactice,  and 
-11  the  unsolv('d  ma  t  iiema  t  i  -  „  I  iroblems  in  tile  fielrl.  rhe  results  i.-resented 
tta'.’  :  rtihablv  save-  efforts  by  otlior  ie-ntists  and  (UKiineers  since  the 
ri  suits  demonstrate  tiic  kind  of  research  tiiaf  '.an  bo  used  in  practice  and 
till-  i  nd  wiiich  is  r'f  mtoslly  i  ht'ort-t  i 'a  I  interest. 


Til.  r'  spo'rodi  i  I  i  ty  for  »  iie  wcr'litsi  and  views  ■  xt  rossed  in  thia  deseriptiv- 
siiinni  ir'.’  I  i'  -  with  Ml-',  .nul  not  witii  tlu  author  f't  tliis  repor*-. 


THEORETICAL  AND  PRACTICAL  ASPECTS  OF  SINGULARITY 

AND  EIGENMODE  EXPANSION  METHODS 

* 

A.  G.  Ramm 


§ 1 .  Int roduction . 

Vast  literati>re  has  been  v;rittcn  on  singularity  and  eiyen- 
mode  expansion  methods  during  the  last  decade.  Engineers  and 
physicists  stimulated  interest  in  the  subject  (see  [1],  [2], 

[3]  and  references  given  in  this  review).  Mathematical  analysis 
of  the  problems  was  initiated  in  [4] ,  [5]  and  was  pushed  con¬ 
siderably  further  by  M.  Agranovich  (see  [1]).  Nevertlieless 
many  questions  in  the  theory  are  open  and  of  considerable  intercr 
to  engineers  and  mathematicians.  The  purpose  of  this  paper  can 
be  summarized  as  follows:  VJe  are  going  to  explain  in  a  simple 
way  the  principal  features  of  the  singularity  and  eigenmodo 
expansion  methods  and  to  formulate  explicitly  1)  what  has  been 
used  by  engineers  with  proof,  2)  what  is  im, port  ant  for  practice 
3)  what  has  been  rigorously  established  and  4)  v;hat  arc  thi- 
un solved  mathematical  problems  in  this  field. 
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The  main  results  obtained  in  scalar  wave  scattering  theory 
were  generalized  to  electromagnetic  wave  scattering  without 
much  difficulty.  That  is  why  and  also  for  siiuplicity  wo  restrict 
ourselves  by  the  presentation  of  the  theory  for  scalar  v/nve 
scattering. 


§2: -  What  are  singularity  and  eigenmode  expansion  methods 
( SEM  and  EEM)  ? 

1.  What  is  EEM? 

Consider  the  problem 


(V^  +  k^)u  =  0  in  U  =  D,  >  0  , 

ulp  =  0  , 

u  =  Uq  +  v  =  exp{ik  (r.)  ,  x)}  +  v 
and  v  satisfies  the  radiation  condition 


(1) 

(2) 

(3) 


?v 


i  k  V  = 


0(1x1 


as 


l-^i 


(4) 


Here  r  is  the  smooth  surface  of  a  finite  obstacle  D. 
If  wo  look  for  a  solution  of  (1)  -  (4)  in  the  form 


u 


cxp{ikr^j,) 


47;r 


f{t)dt,  r 


xt 


xt 


-  X  -  t 


(5) 
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then 


exp(ikr  .  ) 

A(k)f  H  — ^5—: — ^  f(t)  dt  =  -  u.(s)  ,  s  e  r.  (6) 

4.ag^  0 

The  EEM  method  can  now  be  explained  as  follows.  Suppose  that 
the  set  of  eigenvectors  of  operator  A(k) 

A(k)<!)j  =  Aj(k)  ,  j  =  1,2, .  (7) 

2 

forms  a  Riesz  basis  of  L  (r)  =  H.  This  means  that  any  g  e  H 
can  be  expanded  in  the  series 

GO 

j*l  3  3 

and 

<^1  IWII^  i  I  19,1^  i  9  11911^  Cl  >  0  ,  (9) 

^  j=l  J 

2 

where  i|g||  is  the  norm  in  L  (D.  The  inequality  (9)  substi¬ 
tutes  the  Parsoval  equality  for  orthonormal  bases.  A  complete 
system,  in  H  does  not  necessarily  form  a  basis  of  H  (example: 

H  =  L^[0,1],  the  system  =  x^ ,  j  =  0,1,2,...:  Not 

every  g  6  L  [0,1]  can  be  expanded  in  the  series  g{x)  = 

If  the  assumption  made  is  true,  then  equation  (6)  can  be 
solved  by  the  Picard  formula 

CO  u 

f(s)  =  -  y  —  sVCs)  (10) 

3=1  Aj(k)  ^ 

-i- 


where  the  coefficients  are  uniquely  defined  by  the  equality 


u 


0 


I 

3=1 


U  .(h. 

03^3 


(11) 


This  method  of  solution  of  the  scattering  problem  (l)-(4'! 
is  called  EEM.  It  was  used  without  mathematical  analysis 
by  cntjineers  fl],  f3|.  The  ijuostions,  which  immediately 
arise,  can  bo  torinulatod  as  follows:  1)  Is  it.  true  that  the 
nonselfad joint  operator  A(k)  has  eigenvectors  (e.g.  Volterra 
operator  has  no  eigenvectors) ;  2)  Is  it  true  that  the  set  of 
eigenvectors  of  A(k)  forms  a  Riesz  basis  of  H;  3)  suppose 
that  the  set  of  eigenvectors  {h  eigensystem)  of  A(k)  does 
not  form  a  Riesz  basis  of  H.  is  it  true  that  the  root  system 
of  A(k)  forms  a  Riesz  basis  of  H  ? 

Let  us  explain  the  root  system.  Let  A  be  a  linear  operator 
on  H  rtcp  =  X(ti/  (J)  0.  Consider  equation  A(j)j^  ~  1^ 

this  e-’isticn  IS  solvable,  (pj^  is  called  a  root  vector  of  A 
corresponding  to  eigenvalue  A  and  eigenvector  tp  •  If  exists 

consider  equations  A<t)j^  ~  ‘^k-1'*  ^  It  is  known  [7] 

that  only  a  f  i  n  i  to  numtjcr  r  of  root  vectors  , .  .  .  ,  associated 
with  <p  exist.  The  chai  n  ( (ji,  ,),j^ , . .  . ,  tj)^)  is  called  a  Jordan  chain 
with  the  length  r  +  1.  The  union  of  all  root  vectors  of  a 
linear  operator  A  corresponding  to  all  eigenvalues  of  A  is 
called  the  root  system  of  A.  It  is  well  known  from  linear  algebra 
that  the  eigensystem  of  a  nonselfad joint  operator  may  not 
form  a  basis.  For  example  if  the  operator  A  is  an  operator 
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in  IR^  with  the  matrix  A  =  (  )  ,  then  A  has  only  one 

'0  1' 

eigenvector  so  that  the  eigensystem  of  A  does  not  form  a 
2 

basis  of  IR  .  It  is  also  known  that  the  root  system  of  any 
linear  operator  (matrix)  in  Ir”  forms  a  basis  of  IR*'  .  Of 
course  in  IR*'  any  basis  is  Riesz  basis.  In  a  Hilbert  space 
(infinite  dimensional  space)  this  is  not  true.  For  practice 
it  is  important  to  have  af f irmati%’e  answers  for  questions  1) 
and  3) .  Indeed,  if  the  eigensystem  of  A(k)  does  not  form  a 
Riesz  basis  but  its  root  system  forms  a  Riesz  basis  of  H,  then 
it  is  still  possible  to  solve  equation  (6)  using  the  root 
system  of  A. 


2.  What  is  SEM? 

In  order  to  explain  what  SEM  is,  consider  the  problem 

u^t  -  Au  =  0  in  fl  ;  ujj.=  0;  “1^=0  “  “tit=0  ~ 

where  f  6  c”  (a  smooth  function  which  vanishes  for  large  jxj. 

2 

If  G(x,y,-p  )  is  the  Green  function  of  the  problem 

(“A  -»■  p^)G  =  6(x  -  y)  in  fi,  Gjj-  =  0,  Rep  >  0  (13) 

then  the  solution  of  (12)  can  be  written  as 

,c  +  i® 

**  ~  TtFI  I  exp(pt)  U(x,p)  dp,  (14) 


(12) 
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v.’here 


u(x,p) 


2 

G(x,y,-p  )f(y)dy 


(15) 


Suppose  that  u(x,p)  is  meromorphic  in  p  on  the  whole 
con:\plox  plane  p  (this  is  actually  true  [81- [12])  and  suppose 
the  follov-’ing  estimate  is  valid 


u(x,p)  i  £ 


,  a  >  1/2,  c  = const. >  0,  Rep  >  -A, 


1+  IPl 

llmpl  >  , 

v;here  A  >  0  is  arbitrary. 

This  estimate  follov;s  from  the  Lax-Phillips  result  [12] 
and  the  argusucnts  given  in  [9],  [10],  [13].  Under  the  assump¬ 

tions  made,  the  contour  of  integration  in  (14)  can  be  moved  to 
the  left  and  one  gets 


(16) 


N 

u  =  I  c  .t 

j=l  3 


b  -1 

exp(Pjt)  +<?(exp(-Re  Pj^t) ) 


(17) 


where  p^  arc  the  poles  of  G(x,  y,  -p  ) ,  b^  are  the  mul¬ 
tiplicity  of  the  pole  p^.  Expansion  (17)  is  called  the  SEM 
expansion.  Actually  such  expansions  were  known  for  a  long  time 
for  various  concrete  problems  of  matheir.atical  physics  (especially 
in  cases  when  the  solution  can  bo  represented  e.xplicitly  in 
the  form  of  series) .  The  main  difficulty  is  to  prove  estimate 
(1C>)  which  allows  us  to  move  the  contour  on  integration  :n  (14). 
Tf  only  tlio  inormorpliic  nature  of  u(x,p)  as  a  Euncti-'-'  r.f  p 


is  established,  then  no  SE.'l  c:;;.cins ionr  of  the  type  (17)  c£in 
be  proved  in  general  because  there,  is  a  possibility  that  poles 
+  i..',  with  very  snail  e.  >0  and  very  large  ..  can 
exist.  The  general  Mittag-LOff ler  representation  can  not  be 
applied  for  derivation  of  SEM  expansion  (17),  because  this 
representation  uses  special  system  of  expanding  contours,  while 
the  derivation  of  (17)  requires  the  possibility  to  move  our 
particular  contour  (c  -itx.,  c  +ia.)  to  the  left.  From  a  prac¬ 
tical  point  of  view,  the  SEM  expansion  is  used  at  present 
according  to  the  following  scheme: 

Suppose  that  only  n  few  term.s  in  (17)  are  essential,  e.g. 
1-3.  This  will  be  true  if  (Re  p ^ j  »  Rep^  for  j  >  3.  Then 
in  experiments  the  transient  field  u(x,t)  is  measured  and 
each  p  j ,  j  =  1,2,3  is  determined.  It  is  assumed  that  tlse 

2 

location  of  these  complex  F’<^les  of  the  Green  function  G(x,y,-p  ) 
can  give  information  enough  to  identify  the  obstacle  (the 
scatterer  D) .  This  assumption  has  not  been  backed  theoretically 
Nevertheless,  if  there  is  a  finite  set  of  scatters  (say  flying 
targets)  it  is  possible  to  believe  that  a  one  to  one  corre¬ 
spondence  can  be  establJshed  empirically  between  the  scattorers 
and  the  corresponds ng  complex  poles. 

An  intoresi.i  ng  inv^'i  se  prob]  .:..i  can  be  formulated  in  con¬ 
nection  with  thi  ;.  q'.uv:.-c  o);. ; 


Inverse  problem;  Given  a  set  of  co:. itw 


Re  pj  <0,  is  it  j>o. 

Green  functic:  cc:r;  '  s 
{p^^V  Does  fn:.-  set 


yi  c  t('  f  i:;d  u  S'u  t  ; 


win;;  m 


.  .  .  *  O  >■>  ..  I  c  • 


un  iv, ;  i  V  a-. ‘t 


U'hat  restrictions  must  be  imposed  on  the  set  {p  )  ,  Re  o .  <  0 

3  ‘  1 

in  order  that  this  set  will  be  the  sot  of  complex  polos  of  the 
Green  function  of  a  scatterer? 

If  the  scatterer  is  a  star-like  body  and  the  boundax-y  con¬ 
dition  is  the  Dirichlet  condition  then  the  set  fp.}  must 

3 

satisfy  the  condition  |Re  p^j  >  a  In  | Im  |  +  b,  a  >  0  [12]. 
It  seems  that  no  other  information  on  the  problem  is  available. 

From  a  practical  point  of  view  this  problem  may  hot  be  so 
important  as  it  seems.  First,  only  a  fev;  complex  poles  are 
available.  It  seems  hopeless  to  make  any  general  conclusions 
about  the  scaterer  from  this  information  vv’ithout  severe  restric¬ 
tions  on  the  set  of  scatterers.  (For  example,  if  it  is  apriori 
knovm  that  the  scatterer  is  a  ball,  it  is  possible  to  determine 
its  radius  from  the  above  information.)  That  is  why  this  author 
thinks  that  from  a  practical  point  of  view  in  order  to  use  the 
SEM  for  identification  of  scatterers  it  is  more  useful  to  work 
out  tables  of  responses  of  the  typical  scatterers,  then  to  try 
to  develop  a  theory  of  the  posed  (which  is  very  interesting 
from  a  theoretical  point  of  view)  inverse  problem. 


v;ith 


di 


f  f 


3.  The  follov.inr  questions  arise  naturally  in  connection 
the  EEM  a.'in  SEM  r.-'.''thods : 

1)  Docs  the  z'oot  system  ol  the  integral  operators  in 
racticn  theory  (c.g.  opcrator)form  a  Riess  basis  of  K? 

2)  r:'  •  :  p  c  ‘  syrter  f'oincidc-  v;ith  thv'  ciqcn- 


3)  Do  the  complex  poles  of  the  Green  function  depend 

continuously  on  the  obstacle?  In  more  detail:  suppose  that 
Xj  =  ,  1  £  _<  1  are  parametric  equations  of 

r,  Yj  =  +  eZj(t^,t2)/  0  £  0  <  e  <  1 

are  parametric  equations  of  the  surface  of  perturbed  scatterer. 

2 

Let  us  assume  that  Xj(t),  z  ^  (t)  e  C  (A),  t  =  tj^,t2f  A  = 

{t^  _<  ^2/^2  £  !}•  Let  us  fix  an  arbitrary  number  R  >  0  and 
let  p j ,  1  _<  j  _<  r(R)  be  the  complex  poles  of  the  unperturbed 
Green  function  which  lie  in  the  circle  |pj|  <.  R*  Let  Pj{e) 
be  the  complex  poles  of  the  perturbed  Green  function.  Our 
question  can  now  be  formulated  as  follows:  is  it  true  that 
P(e)  Pj»  as  e  0,  uniformly  in  1  £  j  £  n(R)  provided 
that  the  numeration  of  Pj(f;)  is  properly  done? 

4)  How  can  one  calculate  the  complex  poles? 

5)  What  are  sufficient  conditions  for  the  validity  of  SEM 
expansion  (17)? 

6)  Is  it  possible  to  calculate  complex  poles  via  calcula¬ 
tion  of  zeros  of  some  functions? 


§3.  What  has  been  rigorously  established  in  EEM  and  SEM  methods? 

In  this  section  we  give  answers  to  questions  l)-6)  of 
Section  2.3.  No  proofs  will  be  given  but  the  results  obtained 
will  be  formulated  and  references  will  be  given.  Proofs  are 


-9- 


omi  Lted  foe  throf--  rc'ason:;:  ])  they  ai'c  diJ-f.icnilt  fo:;  e!i9  iu'.-i:  i;s , 

2)  they  are  long,  and  3)  they  can  found  in  the  cited  papers. 

1.  In  order  to  formulate  the  answer  to  que.stion  1)  of  Section 

2.3  VIC  must  explain  what  n  Riesz  bas.-' s  v.’ith  brackets  is.  Let 

(f.)  bo  an  orthonorrcal  basis  of  H,  fh  .  }  be  a  co:.,plete  and 
3  3  ■ 

minimal  system  in  H.  (A  complete  system  {h^}  is  called  minimal 
if  the  system  is  not  complete  for  any  k,  k  -  1,2,3... 

In  other  v.’ords  if  v;o  ro;r.ove  any  element  hj^  of  our  system,  we 


obtain  an  incomplet'--  system)  .  Lot  m^  <  rn^  < 


be  an  arbi¬ 


trary  increasing  sequence  of  integers;  F\  is  the  linear  space 

V7ith  the  basis  {f  ,f  .....f  !:  h.  is  the  linear 

’^^‘i-1  'i  -l''  l  ra .  -1  '  3 

space  with  the  basis  {h^^  ,...,h  •  Suppose  that  therc^ 

e.xi.sts  a  linear  bounded  operator  B  v;ith  bounded  B  ’  ,  defined 
on  all  II,  such  that  Bib  -  F  ^ ,  j  -  1,2,...  .  Tlien  L'liO  syscem 
{h^l  is  called  a  Riesz  basis  of  H  with  brac);cts.  This  defi¬ 
nition  i.s  equivalent  (see  [141)  to  the  follov.ring.  I/et  be 

projectors  in  11  onto  lb.  Suppose  th.it  for  any  f  €  II, 


Cl  II  f' 


f  •  2 


1  l  IU’t^II  5  c  ||fj!  ,  c  > 

i=i  ^  ^ 


0  , 


then  the  system  [Im]  is  called  a  Riesz  basis  of  F,  witli 
briickets . 

It  is  proved  thaL  the  root  system  of  op'rntor  h(k)  (see 
formula  (6))  fo’  ms  a  Rio:-./  bt'isi.s  wi  I  h  brackets  [fjj.  The  s-ame 
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is  true  for  the  operator  arising  in  the  exterior  Neumann 
boundary  value  problem  [1],  16].  The  same  is  also  true  for  the 
electrodynamics  scattering  problem  [1]. 

2.  If  the  integral  operator  of  the  diffraction  problem  is 
normal  then  its  eigen  system  coincides  vx'ith  its  root  system. 

This  is  the  case  in  the  problem  (l)-(4)  for  a  spherical  surface 
r  and  for  a  linear  antenna.  First  this  was  observed  in  [5]. 

An  operator  is  called  normal  if  A*A  =  AA*,  where  A*  is  the 
adjoint  operator.  The  condition  AA*  =  A*A  is  the  condition  on 

r  provided  that  the  kernel  of  the  integral  operator  A  is  given, 
(see  [5]  for  details). 

3.  The  answer  to  question  3  from  section  2.3  is  affirmative 
(see  (61  for  details). 

4.  A  general  method  (with  the  proof  of  its  convergence)  for 
calculating  the  complex  poles  of  the  Green's  functions  in  dif¬ 
fraction  and  potential  scattering  theory  was  given  in  (4J,  [15] 
(also  see  [6]).  The  method  can  be  explained  for  the  problem 
with  the  impedance  boundary  condition 

(9^  +  k^)u  =0  in  n  ,  I T  *  ® 

where  h  =  const.#  Re  h  >  0,  N  is  the  outer  normal  to  r  • 

We  look  for  a  solution  of  the  problem  without  sburces  in  the 
form 


-11- 


From  (19)  and  (18)  it  follov.’s  that 


a  =  Qo  , 


whe  re 


(19) 


(20) 


Qo  = 


exp(ikr  )  ,  exp(ikr  ) 

^  o(s)ds  -  h  - a(s)ds  (21) 


Let  {4>j}  be  a  Riesz  basis  of  H  =  L  (F), 


N 


Cj  . 


(22) 


Substituting  (22)  in  (20)  and  multiplying  by  in  H  v;e 

obtain ; 


N 


y  b..(k)c.  =0,  l<i<N;  b..(k)=6-- 

^j=l  11  1  -  -  11  11 


(k)  =  6,^  -  <?;)  .  (23) 


1  1 


This  system  has  nontrivial  solutions  if  and  only  if 


det  b^^ (k)  =  0.  (24) 

The  left  -hand  side  of  (24)  is  an  entire  function  of  k.  Let 

be  its  roots.  Then  there  exist  the  limits  lim  k  =  k 

ra  ..ram 

N-x^ 

and  k^  are  the  polos  of  the  Green  function  corresponding  to 
the  problem  (18).  Moreover,  all  the  complex  poles  can  be  obtained 
by  this  method.  Proofs  are  given  in  [4],  [15],  [6].  From  a 
practicril  jjOLiit  of  viov;  there  are  txvo  nontrivial  points  in  p'or- 
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forming  this  method.  1)  calculation  of  b^j(k)  by  fori.;ulc: 
(23)  and  2)  nuinei'ical  solution  of  equation  (2-1).  For  both 
stops  there  are  methods  available  in  the  literature  on 
numerical  analysis. 

5.  Sufficient  conditions  for  the  validity  of  SEM  expansion 
(17)  v/ere  given  in  Section  2. 

6.  The  set  of  the  complex  polos  of  the  Green  function  of 

the  problem  (1)  -  (4)  coincides  v/ith  the  set  of  the  complex 
zeros  of  the  eigenvalues  operator  A(k)  : 


(25) 


Indeed,  let  G  =  +  ...,  i.e.  z  is  a  pole  of  the  Green 

(k-z) ^ 

function  G,  G|p  =  0, 


G  =  g 
.  -  9G 


-  q(x,t, 


Dm  (t  ,y,k)  dt,  g  = 


exp(ikr^y) 


4Trr 


xy 


aN, 


Multiplying  (26)  by  (k  -  z)  and  taking  k  ■>■  z,  we  get 


(26) 


g{x,t,z)  dt  =  0, 


X  e  r 


(27) 


The  kernel  R(t,y)  is  degenerate, 
such  that 


Thus  a  function  |  0  exists 


g(x,t;z)  (>(t)dt  =  0 


(28) 
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some  of  the  i'liiiet  ioiis 


This  means  that  k  z  is  a  zero  of 

Converse3.y,  if  0  }  0  is  a  solution  of:  {2fi)  ,  thou 

u  =  [  g(>:,t;z)  ^‘;{t)dt  (29) 

is  a  solution  of  the  problem 

(^’^  +  z^)u  =  0  in  o  ,  u  I  .,  =  0  (3  0) 

v/ith  the  outgoing  asynipt.otic  at  infinity.  fJonce  u  =  0  in  Q 

2 

if  7.  is  not  a  pole  of  G.  Since  z  '  is  complex,  u(x)  e  0 
in  D  (as  a  solution  of  the  homogeneous  intei'ior  problem)  . 

By  the  boundary  value  jump  relation  e  0.  This  contradiction 
proves  that  z  is  a  pole  of  G.  A  variational  retliod  for  cal¬ 
culation  of  eigenvalues  of  nonselfacl joi n t  compact  operators 
is  given  in  [6]. 

i 4 .  Open  problems . 

1)  The  inverse  problem  f on-.u’l  ated  in  Section  2  is  of 
interest.  It  is  very  intorcsting  to  have  partial  ansv/ers: 
what  information  about  the  geometry  of  a  scatterer  can  be 
obtained  from  the  locijtion  of  the  complex  poles. 

2)  There  is  a  conjecture  13]  that  tue  com.'lex  poles  of 
the  Green  function  of  th<!  problem  (i)-{4)  for  a  convex  smooth 
compact  boundary  arc  .simple.  It  would  be  interesting  to  [jrove 
it  or  to  give  a  counterexample. 
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3)  It  \.’ould  bo  in  tores  ting  to  perfoiTii  nuinorically  iTi'.tliod 
ciGscribod  in  Section  3.6  in  some  practical  problems. 

4)  In  (IGJ  some!  properties  of  the  purely  real  polos  Ke  p^  <  0, 

Irn  p^  =  0  were  established.  It  w’ould  be  interesting  to  tell 

v.'hat  informa t if;in  about  the  geometry  of  an  obstacle  can  be 

obtained  from  the  location  of  the  purely  real  poles.  Tr,  the  phj'sics 

literature  the  complex  plane  k  =  ip  is  usually  used.  On  this 

plane  the  puredy  real  complex  poles  are  purely  imaginary, 

Re  k  .  =  0 ,  Ira  1:  .  <  0 . 
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